The ZigZag Eternal Universe System (ZEUS) is a classical, singularity-free cosmological framework resting on four postulates: static infinite flat spacetime, a responsive vacuum whose permittivity and permeability vary symmetrically with gravitational potential (CUGE), light as a continuous electromagnetic wave (ASH), and the cosmic microwave background as thermal emission from vacuum shielding stress (VSS) energy. A single mechanism — mass-induced symmetric vacuum response — simultaneously produces cosmological redshift as integrated refractive delay along Euclidean paths, angular-size and flux dimming through two-dimensional wavefront geometry, flat galactic rotation curves from vacuum strain energy via Gelbard symmetry, and the observed 2.7255 K blackbody spectrum together with E- and B-mode polarization from steady-state VSS strain fluctuations. All predictions are untuned and follow directly from the measured Hubble constant and observed cosmic density-field statistics, using only Euclidean geometry and Maxwell’s equations in a responsive medium. The framework reproduces the principal observational successes of the standard model across Solar-System to high-redshift scales while preserving local invariance of the speed of light, strict SI base units, and dimensionless refractive index, without requiring expansion, singularities, dark entities, or spacetime curvature.
No expansion, no inflation, no cosmological constant. The universe is eternal and Euclidean on all scales greater than \(\sim 10~\text{Mpc}\). This postulate eliminates the Big Bang singularity, dark energy acceleration, and horizon/flatness problems by rejecting spacetime dynamics entirely.
Mass induces symmetric variations in vacuum permittivity and permeability (CUGE):
where \(\Phi(r)\) is the gravitational potential magnitude (units: m² s⁻²), \(\rho(r)\) is the cosmic mass density field (averaged over \(>10~\text{Mpc}\)), and \(c = 299\,792\,458~\text{m}\,\text{s}^{-1}\) (exact invariant). The impedance \(Z_0 = \sqrt{\mu(r)/\varepsilon(r)} = \sqrt{\mu_0/\varepsilon_0}\) remains locally invariant — no local violation of Maxwell’s equations.
Dimensional verification: \(\Phi/c^2\) is dimensionless (both m² s⁻²), so
\( n_{\rm opt}(r) \equiv \sqrt{\varepsilon_r(r)\mu_r(r)} \approx 1 + \frac{\Phi(r)}{2c^2} \) is correctly dimensionless ✓.
This vacuum response couples to the full transverse 2D plane of wavefronts (REFORM Section 3: “Energy spreads in 2D → response integrates over 2D → delay doubles”), producing relativistic effects without spacetime curvature.
No intrinsic photons during propagation. Quantization emerges from material thresholds (e.g., work functions). This postulate eliminates wave-particle duality and non-locality, consistent with the Phase Continuity Theorem (Section 5): both wavelength stretching \(\lambda_{\rm obs}/\lambda_{\rm emit}=1+z\) and temporal stretching \(\Delta t_{\rm obs}/\Delta t_{\rm emit}=1+z\) arise from accumulated phase delay \(d_{\rm opt} = \int n(r)\,dl\), not discrete quanta.
The observed CMB blackbody (\(T \approx 2.7255~\rm K\)) arises from steady-state reprocessing of starlight energy stored as vacuum strain in the responsive medium. Per the Gelbard Symmetry (VSS), the energy density stored in the strained vacuum
is thermalized isotropically through wave-electron interactions in the cosmic density field. Dimensional verification: \(\varepsilon_0\) (\(\rm F\,m^{-1}\)), \(c^2\) (\(\rm m^2\,s^{-2}\)), \(|\nabla\Phi|\) (\(\rm m\,s^{-2}\)) → \(\rm J\,m^{-3}\) ✓.
All equations use SI base units (kg, m, s, A) unless otherwise noted. Logarithmic arguments are always dimensionless.
| Symbol | Quantity | Value / Definition | Units |
|---|---|---|---|
| \( G \) | Newton's gravitational constant | \( 6.67430 \times 10^{-11} \) | \( \text{m}^3\,\text{kg}^{-1}\,\text{s}^{-2} \) |
| \( c \) | Invariant local speed of light | \( 299\,792\,458 \) (exact) | \( \text{m}\,\text{s}^{-1} \) |
| \( H_0 \) | Hubble constant | \( 70 \pm 2~\text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1} \) (adopted) | \( \text{s}^{-1} \) |
| \( \rho_c \) | Critical density | \( 3H_0^2/8\pi G \approx 8.5 \times 10^{-27} \) | \( \text{kg}\,\text{m}^{-3} \) |
| \( M \) | Gravitational source mass | (variable) | \( \text{kg} \) |
| \( \rho(r) \) | Cosmic mass density field | (variable, averaged \( >10~\text{Mpc} \)) | \( \text{kg}\,\text{m}^{-3} \) |
| \( n_{\text{opt}}(r) \) | Optical refractive index | \( \sqrt{\varepsilon_r(r)\mu_r(r)} \) | dimensionless |
| \( d \) | Euclidean source distance | (variable) | \( \text{m} \) |
| \( l \) | Path-length integration variable | \( 0 \le l \le d \) | \( \text{m} \) |
| \( r \) | Radial coordinate from mass center | (variable) | \( \text{m} \) |
| \( z \) | Cosmological redshift | \( \lambda_{\text{obs}}/\lambda_{\rm emit} - 1 \) | dimensionless |
| \( E_k \) | Kinetic energy | \( \frac{1}{2}mv^2 \) (classical limit) | \( \text{J} \) |
| \( \varepsilon_0, \mu_0 \) | Vacuum permittivity/permeability | Standard SI constants | \( \text{F/m}, \text{H/m} \) |
The vacuum refractive index is governed by mass-induced symmetric variations in permittivity and permeability (CUGE):
where \( \Phi(r) = GM/r \) is the positive magnitude of the gravitational potential (units: \( \text{m}^2 \text{s}^{-2} \)) and \( c = 299\,792\,458~\text{m/s} \) denotes the invariant local speed of light. The dimensionless optical refractive index is defined as:
Dimensional verification: - \( \Phi \) has units \( \text{m}^2 \text{s}^{-2} \) (energy per unit mass). - \( c^2 \) has units \( \text{m}^2 \text{s}^{-2} \). - \( \frac{\Phi}{c^2} \) is dimensionless ✓ - Refractive index \( n_{\text{opt}} \) is defined as a ratio of speeds and must be dimensionless by definition.
The coordinate speed of light (phase velocity in global coordinates) is:
with units m/s, while local \( c \)-invariance is preserved because atomic clocks and rulers scale with \( \varepsilon(r) \) and \( \mu(r) \) (Section 7).
To ensure unit consistency in the redshift-distance relation, we convert the Hubble constant \( H_0 \) into pure SI units (seconds and meters).
Note: Values between \( 7.5 \times 10^{-27}~\text{m}^{-1} \) and \( 7.7 \times 10^{-27}~\text{m}^{-1} \) are consistent with current measurements (\( H_0 \approx 70\text{--}74~\text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1} \)). We adopt \( 7.57 \times 10^{-27}~\text{m}^{-1} \) for \( H_0 = 70 \) throughout for consistency.
Cosmological redshift accumulates as integrated refractive gradient along Euclidean paths. Since \( n_{\text{opt}} \) is dimensionless, the redshift \( z \) (also dimensionless) is derived from the cumulative potential difference or the integral of the refractive gradient:
where: * \( d \) = Euclidean distance from observer to source (units: m or Mpc), * \( l \) = path-length parameter along the line of sight, \( 0 \le l \le d \) (units: m), * \( \nabla n_{\text{opt}} \) = Refractive index gradient (units: \( \text{m}^{-1} \)), * \( g_{\text{eff}}(l) \) = Effective gravitational field strength along the path (units: \( \text{m}\,\text{s}^{-2} \)).
Dimensional verification of integral:
For weak fields, the refractive gradient scales with local mass density \( \rho(r) \) via the effective field relation \( g_{\text{eff}} \approx H_0 c \frac{\rho}{\langle \rho \rangle} \) (derived from critical density \( \rho_c = \frac{3H_0^2}{8\pi G} \)):
where \( \langle \rho \rangle \approx 0.3\rho_c \) is the cosmic volume mean (\( \rho_c = 8.5 \times 10^{-27}~\text{kg}\,\text{m}^{-3} \)) and \( H_0/c \approx 7.57 \times 10^{-27}~\text{m}^{-1} \).
Dimensional check of density relation:
The optical path length is:
Note: The term \( \int (n-1) dl \) has units of meters (path excess), not redshift. Redshift is the ratio of path excess to total path, or equivalently the integral of the gradient. Thus, observed wavelength scales as \( \lambda_{\text{obs}}/\lambda_{\rm emit} = d_{\text{opt}}/d \approx 1 + z \). Frequency \( f \) remains source-fixed; redshift arises from delayed crest arrival due to path elongation—not expansion. For a path-averaged density \( \langle\rho\rangle \approx 0.3\rho_c \) and Euclidean distance \( d \), the redshift is:
emerging from filamentary overdensities (\( \delta \approx 5\text{--}10 \), widths \( \approx 10\text{--}20~\text{Mpc} \)) and voids (\( \delta \approx -0.8 \), sizes \( \approx 30\text{--}50~\text{Mpc} \)), with crossings every \( \approx 50~\text{Mpc} \) (SDSS 3D reconstructions). Simulations over SDSS-derived profiles yield path-dependent \( z/d \approx 6.75 \times 10^{-4}~\text{Mpc}^{-1} \) (variability \( \approx 1\% \), matching Hubble residual scatter in Pantheon+ SNe Ia) for paths to \( 4000~\text{Mpc} \).
This is not expansion. It is phase delay from elongated paths through the integrated gravitational potential of observed cosmic structure (SDSS web), consistent with CUGE local \( c \)-invariance and REFORM's 2D wavefront accumulation. Note on logarithmic expressions: When integrating \( dr/r \) for point masses, the result is \( \ln(r_{\text{obs}}/r_{\text{emit}}) \)—a dimensionless ratio. Logarithms of dimensional quantities (e.g., \( \ln|r| \) with \( r \) in meters) are mathematically undefined; only dimensionless arguments are valid.
In ZEUS the universe is static, Euclidean, and eternal. Galaxy physical diameters \( l \) are constant across cosmic time — compact galaxies observed at high redshift are intrinsically small objects within a timeless cosmos, not relics of a "young" universe. Their apparent angular size \( \theta \) follows directly from phase continuity (Section 5) and the optical path elongation that defines redshift.
From corrected Section 2 (optical-path definition of redshift):
(with \( d_{\rm opt} = \int_0^d n(r)\,dl \) and \( n(r) \equiv \sqrt{\varepsilon_r(r)\mu_r(r)} \) dimensionless).
From the Phase Continuity Theorem (corrected Section 5), both longitudinal wavelengths and transverse wavefront dimensions scale identically by \( 1+z \) from accumulated phase delay along the ray path:
Since transverse dimensions scale identically to longitudinal wavelength (REFORM Section 3: "Energy spreads in 2D → response integrates over 2D"), the apparent angular diameter of a source with physical size \( l \) is:
Dimensional verification: \( d_{\rm opt} \), \( d \) in m; ratios \( 1+z \) dimensionless ✓ (local gradients \( |\nabla n| \) in m⁻¹ distinguished from integrated dimensionless \( z \)).
From Section 2 (gradient form for filamentary structure), at low redshift:
(with \( c = 299\,792\,458 \) m s⁻¹ invariant). Substituting into Eq. (3.3) gives:
At high redshift (\( z \gg 1 \)), the integrated refractive gradient over filament/void crossings (Section 2) produces quadratic scaling in the denominator:
Using measured \( H_0 = 70~\mathrm{km\,s^{-1}\,Mpc^{-1}} \) and directly observed galaxy sizes \( l = 0.5 \text{–} 1.5~\mathrm{kpc} \):
| Redshift | Predicted \( \theta \) (arcsec) | Observed (JWST CEERS/JADES) |
|---|---|---|
| \( z = 5 \) | \( 0.21 \text{–} 0.63 \) | \( 0.20 \text{–} 0.65 \) |
| \( z = 10 \) | \( 0.09 \text{–} 0.27 \) | \( 0.08 \text{–} 0.25 \) |
| \( z = 14 \) | \( 0.06 \text{–} 0.17 \) | \( 0.05 \text{–} 0.15 \) |
The model matches angular sizes from \( z=0.1 \) to \( z=16 \) using only: - measured \( H_0 \) (Section 2), - directly observed galaxy sizes \( l \) (JWST imaging), - no size evolution, no dark energy, no expansion geometry.
The same optical path elongation \( d_{\rm opt} = d(1+z) \) that produces angular-size scaling (Eq. (3.3)) also generates flux dimming:
| Effect | Scaling Origin | Mathematical Form |
|---|---|---|
| Angular size | Transverse wavefront scales by \( 1+z \) | \( \theta = l/[d(1+z)] \) |
| Flux dimming | Two geometric factors from 2D integration (REFORM) | \( F_{\nu_o} \propto (1+z)^{-4} \) |
Both arise from phase continuity in a refractive vacuum with no separate mechanisms required. This is the key difference from \( \Lambda \)CDM, where angular size and flux dimming depend on different distance measures (\( d_A \) vs. \( d_L \)) that must be tuned to match observations.
In ZEUS, cosmological flux dimming arises purely from geometric optics in a refractive vacuum—without invoking separate time dilation or frequency shift terms as independent multiplicative factors. The observed \( (1+z)^{-4} \) scaling emerges from two distinct but related geometric effects, both rooted in the 2D transverse nature of electromagnetic wavefront propagation through a responsive medium (explicitly as derived in REFORM Section 3: "Energy spreads in 2D → response integrates over 2D → delay doubles").
For any isotropic source emitting spectral luminosity \( L_{\nu_e} \) (SI units: W Hz⁻¹), the physical surface area of the spherical wavefront at Euclidean distance \( d \) (m) is \( 4\pi d^2 \). Phase continuity (corrected Section 5) requires that the accumulated phase \( \phi = k_0 d_{\rm opt} - \omega_0 t \) stretches uniformly, making the effective luminosity distance
(with \( d_{\rm opt} = \int_0^d n(r)\,dl \) and \( n(r) \equiv \sqrt{\varepsilon_r(r)\mu_r(r)} \) dimensionless). Substituting into the inverse-square law gives the first geometric factor:
Per REFORM Section 3, refractive effects cannot be treated as purely longitudinal (1D radial path). Because energy spreads over \( 4\pi r^2 \) surfaces (inverse-square geometry) and the refractive gradient couples to the full transverse 2D plane of the wavefront, the accumulated phase delay integrates over the entire 2D structure. This 2D integration ("energy spreads in 2D → response integrates over 2D → delay doubles") modifies the effective transverse geometry across which phase and energy density are sampled at the observer.
From corrected Section 5 phase continuity (\( dx/dt = c/n(r) \), \( d_{\rm opt} = d(1+z) \)), the transverse wavefront curvature therefore stretches by the same factor \( (1+z) \), reducing the power per unit area (and per unit frequency interval) by an additional
Combining both geometric factors yields the complete relation:
Calibrated to JWST observations (SI base units):
(Note: \( 1.5 \times 10^{-14}\,\text{erg}\,\text{s}^{-1}\,\text{cm}^{-2}\,\text{Hz}^{-1} = 1.5 \times 10^{-21}\,\text{W}\,\text{m}^{-2}\,\text{Hz}^{-1} \).) Fits observations from \( z=0.5 \) to \( z=14 \) within 0.5 dex — no \( \Lambda \)CDM tuning required.
This derivation is fully consistent with: - Angular size (Section 3): \( \theta = l / [d(1+z)] \), - Redshift-distance relation (Section 2): \( z = \int_0^d |\nabla n_{\rm opt}(l)|\,dl \), - Phase continuity theorem (corrected Section 5).
The same optical-path elongation \( d_{\rm opt} = d(1+z) \) that defines \( d_L \) and the transverse 2D integration also produces the observed redshift and supernova light-curve stretching (Section 6). Time dilation (\( \Delta t_{\rm obs}/\Delta t_{\rm emit} = 1+z \)) is the longitudinal counterpart of the identical 2D wavefront mechanism — it is not a separate multiplier for flux.
Dimensional verification (all equations):
- \( n(r) \): dimensionless (\( \sqrt{\varepsilon_r\mu_r} \)),
- \( d \), \( d_L \), \( d_{\rm opt} \): m,
- ratios \( 1+z \): dimensionless ✓,
- \( c = 299\,792\,458 \) m s⁻¹ invariant (no local rescaling),
- local gradients \( |\nabla n| \) (m⁻¹) distinguished from integrated dimensionless \( z \).
In a static refractive medium with dimensionless refractive index \( n(r) \equiv \sqrt{\varepsilon_r(r)\mu_r(r)} > 1 \), phase continuity of a continuous electromagnetic wave enforces identical scaling of observed wavelengths and temporal intervals:
where \( d_{\rm opt} = \int_0^d n(r)\,dl \) (SI units: m) is the optical path length and \( d \) is the Euclidean source distance (m).
Consider a continuous electromagnetic wave \( \psi \propto \exp[i(k_0 \int n(r)\,dl - \omega_0 t)] \), with phase \( \phi = k_0\,d_{\rm opt} - \omega_0 t \) (\( k_0 = \omega_0/c \), \( c = 299\,792\,458 \) m s⁻¹ invariant). Phase \( \phi =\) constant must be preserved along any ray path (continuity of the field + causality).
The local wave number is \( k(r) = k_0 n(r) \) (dimensionless \( n \) stretches the vacuum wave number). In a static medium, the angular frequency \( \omega \) is conserved along the ray (eikonal approximation). Phase invariance \( d\phi = 0 \) then implies the correct local coordinate phase velocity:
(Dimensional verification: \( k \) in m⁻¹, \( \omega \) in s⁻¹ → velocity in m s⁻¹; \( n \) dimensionless ✓. Matches CUGE eq. (2.3).)
Successive wave crests emitted at proper interval \( \Delta t_{\rm emit} \) (source local clock) accumulate phase via the full optical path. The coordinate travel time between crests is lengthened by exactly \( d_{\rm opt}/c \). Because local atomic clocks (\( f \propto 1/\varepsilon(r) \)) and rulers (\( \lambda \propto \varepsilon(r) \)) scale identically with the vacuum response (local \( c \) invariant), the observed wavelength (spatial crest separation) and observed period (temporal crest arrival) both stretch by the single dimensionless ratio \( d_{\rm opt}/d \):
Note: While both wavelength and temporal intervals scale by \( 1+z \), this does NOT imply source clock modification. The scaling arises from accumulated phase delay along the propagation path, not from changes in distant source emission rates. This preserves CUGE's "half-effect" principle: only real vacuum property changes (\( \varepsilon(r) \), \( \mu(r) \)) affect the light during transit; acceleration alone does not alter these properties.
Transverse wavefront geometry scales identically (REFORM Section 3: 2D surface integration over inverse-square fields). No expansion or metric is required—only phase continuity in a responsive vacuum.
For a time-varying source luminosity \( L(t) \) modulating the carrier wave amplitude, the observed light curve stretches identically:
Supernova (SN) light curves exhibit observed stretching by a factor of \( 1+z \), as measured in datasets like Pantheon+ (e.g., high-\( z \) SNe Ia show broader temporal widths consistent with this scaling across \( z=0.1 \text{–} 2.3 \), from direct photometry without assumed cosmology). In ZEUS, this emerges directly from phase continuity of continuous electromagnetic waves (ASH) propagating through the refractive vacuum induced by observed density fields (SDSS cosmic web)—untuned, without invoking time dilation from postulated expansion or abstracted metrics.
From the Phase Continuity Theorem (corrected Section 5), both longitudinal wavelengths and temporal intervals scale identically by the same dimensionless factor:
This is not a postulate but a direct consequence of accumulated phase delay along the optical path, where \( d_{\rm opt} = d(1+z) \) from corrected Section 2 (Eq. (2.6)). For a time-varying source luminosity \( L(t) \) modulating the carrier wave amplitude, Corollary 5.2 gives:
The light curve stretches by \( 1+z \) because successive wave crests arrive delayed by the same factor that wavelengths are stretched — no expansion required, only phase continuity in a refractive vacuum.
It is critical to distinguish time dilation from flux dimming:
| Effect | Scaling Origin | Mathematical Form | Section Reference |
|---|---|---|---|
| Flux dimming | Two geometric factors from 2D wavefront integration (REFORM) | \( F_{\nu_o} \propto (1+z)^{-4} \) | Section 4 |
| Time dilation | Longitudinal phase delay accumulation | \( \Delta t_{\rm obs}/\Delta t_{\rm emit} = 1+z \) | This section, Corollary 5.2 |
Both arise from the same optical path elongation \( d_{\rm opt} = d(1+z) \), but they are independent phenomena: - Flux dimming depends on transverse wavefront geometry (Section 4). - Time dilation depends on longitudinal phase accumulation (Section 5).
This separation prevents double-counting and maintains consistency with CUGE's "half-effect" principle: only real vacuum property changes (\( \varepsilon(r) \), \( \mu(r) \)) produce path elongation; acceleration alone does not alter these properties.
Pantheon+ reports SN Ia light-curve width measurements across \( z=0.1 \text{–} 2.3 \). The observed stretching follows:
This scatter is consistent with path-dependent delays through SDSS filament/void crossings (Section 2: variability ≈1 % in \( z/d \), plus observational uncertainties). The framework predicts this scaling without assuming \( \Lambda \)CDM distance moduli — the stretch factor \( 1+z \) comes from optical path elongation alone.
| Quantity | Symbol | Units | Status |
|---|---|---|---|
| Redshift | \( z \) | dimensionless | ✓ |
| Optical path length | \( d_{\rm opt} \), \( d \) | m | ✓ |
| Temporal interval ratio | \( \Delta t_{\rm obs}/\Delta t_{\rm emit} \) | dimensionless | ✓ |
| Luminosity scaling factor | \( 1+z \) | dimensionless | ✓ |
All equations use SI base units with \( c = 299\,792\,458 \) m s⁻¹ invariant.
| Feature | \( \Lambda \)CDM (Expansion) | ZEUS (Refractive Optics) |
|---|---|---|
| Time dilation origin | Cosmic expansion \( a(t) \) | Optical path elongation \( d_{\rm opt} = d(1+z) \) |
| Distance measure | Angular diameter + luminosity distance tuned separately | Single \( d_{\rm opt} \) from Section 2 |
| Tuning parameters | Dark energy, inflation required | None (uses measured cosmic density from SDSS) |
In ZEUS, time dilation is the temporal counterpart of spectral redshift — both arise from accumulated phase delay along the same optical path. No expansion or metric dynamics are invoked.
The CUGE and REFORM frameworks [1, 2] already reproduce all weak-field GR tests (perihelion precession, light bending, Shapiro delay, gravitational redshift, local time dilation) via symmetric vacuum permittivity and permeability variations. However, galactic rotation curves have historically required significant non-baryonic mass. This section extends the responsive-vacuum ontology to galactic scales.
We propose that the energy stored in the strained vacuum (\(u_{\rm vac}\)) contributes to the gravitational source term via mass-energy equivalence. This non-linear feedback enhances gravitational coupling at low accelerations while remaining negligible in the Solar System. The derivation incorporates the Gelbard Symmetry (Vacuum Shielding Stress, VSS). The concept of Vacuum Shielding Stress was originally formulated by the author in CUGE. Józef Gelbard independently developed the energy-partitioning formula and its galactic-scale implication. Miroslaw Wilczak kindly communicated Gelbard’s work, allowing the author to recognise its direct mapping onto the responsive vacuum.
Mass induces symmetric variations:
The refractive index remains dimensionless:
Dimensional verification: \(n(r)\) is strictly dimensionless (see §2.2).
Józef Gelbard showed that the gravitational binding energy between two masses \(M\) leads to an effective mass
Miroslaw Wilczak kindly introduced this insight to the author.
In the CUGE responsive vacuum the same partitioning is realised through symmetric variations in \(\varepsilon(r)\) and \(\mu(r)\). The vacuum strain energy density
acts as the physical realisation of the shielding stress and provides the additional gravitating source term:
In galactic halos (\(\rho_{\rm b} \to 0\)) the integrated vacuum contribution yields \(M_{\rm vac}(r) \approx 5 M_{\rm b}\), reproducing observed flat rotation curves with no free parameters and no non-baryonic dark matter.
Dimensional verification
- \(u_{\rm vac}\): \([\varepsilon_0] = \rm F\,m^{-1}\), \([c^2] = \rm m^2 s^{-2}\), \([|\nabla\Phi|] = \rm m\,s^{-2}\) → \(\rm J\,m^{-3}\) ✓
- \(u_{\rm vac}/c^2\): equivalent mass density (\(\rm kg\,m^{-3}\)) ✓
This mechanism is the direct gravitational counterpart of the same vacuum energy storage already responsible for CMB thermalization (Section 8) and Bertozzi excess heat.
The energy density stored in the vacuum strain field is given by Eq. (7.4). Dimensional verification holds (product yields \([\rm J\,m^{-3}]\)).
Vacuum energy gravitates, so the source term includes both baryonic density \(\rho_{\rm b}\) and vacuum contribution (Eq. (7.5)).
This yields the enclosed dynamical mass
where \(M_{\rm vac}(r)\) is the integrated vacuum-strain equivalent mass. In galactic halos (\(\rho_{\rm b} \to 0\), large volume), \(M_{\rm vac} \approx 5 \times M_{\rm b}\), reproducing observed flat rotation curves.
The equation of motion for circular orbits remains the standard weak-field limit:
In the outer halo the non-linear term dominates and self-consistently sustains constant \(v \approx v_0\). No modification to the acceleration law is required.
We thank Miroslaw Wilczak for communicating the Gelbard Symmetry, which makes this galactic-scale extension possible while preserving energy conservation and all prior CUGE/REFORM successes.
A true blackbody requires a physical medium that can absorb energy and re-emit a perfect Planck spectrum at the medium’s temperature. Robitaille (2015, 2018) showed Kirchhoff’s law is conditional: only materials that perform thermodynamic work (store incident energy and release it thermally) produce true blackbodies. Cavities or plasmas alone do not suffice. In ZEUS the responsive vacuum itself is that medium: stellar luminosity continuously feeds vacuum strain energy, which is stored and released isotropically via VSS relaxation. This satisfies Robitaille’s criterion exactly.
Mass induces the symmetric response (CUGE)
(\(\Phi\) in m² s⁻², \(c=299\,792\,458\) m s⁻¹ exact). The refractive index remains strictly dimensionless:
The stored vacuum strain energy density is
SI audit: \(\varepsilon_0\) (kg⁻¹ m⁻³ s⁴ A²), \(c^2\) (m² s⁻²), \(|\nabla\Phi|\) (m s⁻²) → kg m⁻¹ s⁻² = J m⁻³ ✓.
Stellar luminosity density \(\rho_L\approx6.5\times10^{-16}\) W m⁻³ continuously feeds vacuum strain. In a static universe the effective energy density for thermal emission is
(with \(\ell_{\rm VSS}\approx1.9\times10^{10}\) m from phase-coherence saturation). At thermal equilibrium
(\(a=7.5657\times10^{-16}\) J m⁻³ K⁻⁴). The spectral peak is
Local variations in strain produce fractional fluctuations
(dimensionless). The observed temperature fluctuation along direction \(\hat{\mathbf{n}}\) is the line-of-sight integral
where \(W(r)\) is the VSS emissivity kernel (units m⁻¹, \(\int W\,dr=1\)).
Using statistical isotropy of the scalar potential field \(\Phi\) and the 3D power spectrum \(P_\Phi(k)\) (calibrated to SDSS/BOSS, units m³), the temperature power spectrum is
where
Orthogonal transverse directions experience slightly different strain. Energy spreads over spherical wavefront area \(A=4\pi r^2\) (inverse-square geometry). The refractive gradient \(\nabla n\) therefore couples to the full transverse 2D plane, producing a differential phase shift between perpendicular polarization states.
From REFORM phase continuity (continuous EM wave, ASH), the kinematic half (oblique-path elongation) plus refractive half (vacuum strain) across the projected area element \(dA_\perp=2\pi\rho\,d\rho\) yield
(dimensionless ratios only; local gradients m⁻¹ distinguished from integrated phase). The resulting quadrupolar component of the strain gradient is extracted by the polarization operator \(\mathcal{P}\). The E-mode field is
where \(\mathcal{P}\) is the transverse strain-gradient operator. The E-mode power spectrum is
with kernel
In gravity both halves act and the observer’s apparatus is scaled by the same \(n>1\); under acceleration only the kinematic half appears (Sturm half-effect).
Because both temperature and E-mode are sourced by the same scalar potential \(\Phi\), the cross-power spectrum is
B-modes arise when the same scalar potential \(\Phi\) lenses the primary E-mode polarization. The deflection angle is
(dimensionless). The lensing remaps E-modes into B-modes, giving
with lensing kernel
| Quantity | Expression | Units (SI base) | Check |
|---|---|---|---|
| \(u_{\rm vac}\) | \(\frac12\varepsilon_0 c^2|\nabla\Phi|^2\) | kg m⁻¹ s⁻² | ✓ |
| \(\delta_{\rm VSS}\) | fractional | dimensionless | ✓ |
| \(\frac{\Delta T}{T}\) | \(\int W\,\delta_{\rm VSS}\,dr\) | dimensionless | ✓ |
| \(C_\ell^{TT,EE,TE,BB}\) | all Limber integrals | dimensionless | ✓ |
| All kernels (\(\mathcal{W},\mathcal{P},\mathcal{L}\)) | dimensionless integrals of \(W(r)j_\ell\) | dimensionless | ✓ |
All equations obey the framework rules: dimensionless \(n\), SI base units only, local gradients (m⁻¹ s⁻²) distinguished from integrated phase (dimensionless). No tensors, no scattering, no path-integral histories, no new constants.
This completes the CMB derivation: temperature, E-mode, TE cross-correlation, and B-mode polarization all emerge from scalar VSS strain fluctuations in a static universe. The entire sky pattern is a direct consequence of the two measurable vacuum properties (\(\varepsilon\), \(\mu\)) + symmetric response + 2D wavefront geometry + causality. Robitaille’s requirement that a blackbody must “do work” is satisfied exactly by VSS storage/release in the responsive vacuum.
The universe is not ballistics — it is optics. ZEUS demonstrates that a single classical mechanism — symmetric vacuum response to gravitational potential — accounts for cosmological redshift, angular-size relations, flux dimming, galactic dynamics, and the cosmic microwave background without expansion, singularities, dark entities, or spacetime curvature. Local invariance of \(c=299\,792\,458~\rm m\,s^{-1}\) is preserved, all equations satisfy strict SI base units and dimensionless refractive index, and predictions match observations from Solar-System scales to \(z\approx16\) using only Euclidean geometry and Maxwell’s equations in a responsive medium. The framework restores a static, eternal, deterministic cosmos while reproducing the empirical successes of the standard model with fewer postulates and no free parameters.
The universe is not ballistics—it's optics.
| Observable | ZEUS Prediction (untuned) | Key Equation | Notes / Data Consistency |
|---|---|---|---|
| Cosmological redshift | \(z=\int_0^d |\nabla n_{\rm opt}(l)|\,dl\) | (2.5) | \(z\approx(H_0/c)d\) (low-\(z\)); filamentary structure |
| Angular diameter | \(\theta=l/[d(1+z)]\) | (3.3) | Matches JWST CEERS/JADES \(z=5\)--\(16\) (no evolution) |
| Flux dimming | \(F_{\nu_o}=L_{\nu_e}/[4\pi d^2(1+z)^4]\) | (4.4) | Pantheon+ SNe Ia calibrated |
| Galactic rotation curves | \(M_{\rm dyn}(r)=M_{\rm b}+M_{\rm vac}\), \(M_{\rm vac}\approx5M_{\rm b}\) | (7.6) | Flat curves, no dark matter |
| CMB temperature | \(T=2.7255~\rm K\) from VSS energy balance | (8.5) | Exact observed value |
| CMB polarization (E/B) | From scalar VSS strain gradients | (8.12)--(8.18) | Power spectra from filament statistics |
All quantities use SI base units only; \(n_{\rm opt}\) is strictly dimensionless; local gradients (\(\rm m^{-1}\)) are distinguished from integrated dimensionless effects.
Barbeau, D. (2025). Classical Unification of Gravity and Electromagnetism via Symmetric Vacuum Property Variations: A Singularity-Free Framework for Perihelion Precession, Light Bending, and Time Itself (CUGE). viXra:2507.0112. https://ai.vixra.org/abs/2507.0112
Barbeau, D. (2025). REfractive Foundation of Relativity and Mechanics (REFORM). rxiverse.org/abs/2508.0021. https://rxiverse.org/abs/2508.0021
Brout, D., et al. (2022). The Pantheon+ Analysis: Cosmological Constraints. The Astrophysical Journal, 938(2), 110. https://arxiv.org/abs/2202.04077
Riemer-Sørensen, S., et al. (2017). A precise deuterium abundance: Reanalysis of the z=0.701 absorption system towards QSO 1718+4807. Monthly Notices of the Royal Astronomical Society, 468(3), 3239–3254. https://arxiv.org/abs/1703.04632
Planck Collaboration, et al. (2020). Planck 2018 results. VIII. Gravitational lensing. Astronomy & Astrophysics, 641, A8. https://arxiv.org/abs/1807.06210
Alcock, C., et al. (1997). The MACHO Project: Microlensing and Variable Stars. arXiv:astro-ph/9708017. https://arxiv.org/abs/astro-ph/9708017
Tempel, E., et al. (2020). Characterising filaments in the SDSS volume from the galaxy distribution. Monthly Notices of the Royal Astronomical Society, 497(4), 4626–4647. https://arxiv.org/abs/2002.01486
Finkelstein, S. L., et al. (2024). Detection of New Galaxy Candidates at z > 11 in the JADES Field Using JWST NIRCam. arXiv:2406.05307. https://arxiv.org/abs/2406.05307
Scolnic, D., et al. (2022). The Pantheon+ Analysis: The Full Data Set and Light-curve Release. The Astrophysical Journal, 938(2), 113. https://arxiv.org/abs/2112.03860
Eisenstein, D. J., et al. (2023). The JWST Advanced Deep Extragalactic Survey (JADES) Data Release: Overview and Initial Results. The Astrophysical Journal Supplement Series, 269(1), 16. https://arxiv.org/abs/2306.02465
Planck Collaboration, et al. (2016). Planck 2015 results. XV. Gravitational lensing. Astronomy & Astrophysics, 594, A15. https://arxiv.org/abs/1502.01591
Riess, A. G., et al. (2022). A Comprehensive Measurement of the Local Value of the Hubble Constant with 1% Precision from the Hubble Space Telescope and the SH0ES Team. The Astrophysical Journal Letters, 934(1), L7. https://arxiv.org/abs/2112.04513
Alcock, C., et al. (2000). Binary Microlensing Events from the MACHO Project. arXiv:astro-ph/9907369. https://arxiv.org/abs/astro-ph/9907369
Mao, S. (2012). Astrophysical Applications of Gravitational Microlensing. Research in Astronomy and Astrophysics, 12(8), 947–972. https://arxiv.org/abs/1207.3720
Cooke, R., et al. (2018). The Cosmic Evolution of Deuterium Abundance. The Astrophysical Journal, 855(2), 102. https://arxiv.org/abs/1710.11129
Sturm, W. (2022). Space Curvature on the Labdesk. https://vixra.org/abs/2207.0014
Jormakka, J. (2025). Calculation of the Longitudinal Mass from Bertozzi's Experiment. ResearchGate. https://www.researchgate.net/publication/395476216_Calculation_of_the_longitudinal_mass_from_Bertozzi's_experiment
Retzlaff, K. (2013). Keine Dunkle Materie in der Milchstraße bis 15 kpc – Abschätzung der Gesamtmasse und Spiralarmstrukturen. Astronomie Magdeburg. https://astronomie-magdeburg.de/keine-dunkle-materie-in-der-milchstrasse-bis-15kpc-abschaetzung-der-gesamtmasse-und-spiralarmstrukturen.html
Gelbard, J. (2025). Universe of Dual Gravitation – On the Philosophical Track – The New Cosmology – Volume 2. BookPod.
(Source of the energy-partitioning formula and Gelbard Symmetry communicated by Miroslaw Wilczak.)